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Luminescence as a mechanism for terahertz emission from femtosecond laser-induced plasmas is 
studied. By using a fully microscopic theory, Coulomb scattering between electrons and ions is 
shown to lead to luminescence even for a spatially homogeneous plasma. The spectral features 
introduced by the rod geometry of laser-induced plasma channels in air are discussed on the basis 
of a generalized mode-function analysis. 
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The past decade has seen increasing activity in the gen- 
eration and detection of far-infrared or terahertz (THz) 
radiation with accompanying increase in applications, in- 
cludingtime-domain THz spectroscopy, and THz tomog- 
raphy [ij. Current THz sources have limited spatial ex- 
tent (/^m to mm), such as transient antenna emitters, 
sources based on optical rectification, and quantum cas- 
cade lasers. A notable exception is the observation of 
THz radiation emitted from plasma channels induced in 



air by femtosecond (fs) infrared filaments ^ 3. 4J , or light 
strings, that can extend over centimeters HiBlMSl- The 
extension of these results to other media offers the pos- 
sibility of a new variety of extended THz sources. 

In addition to the technological importance of light 
strings for realizing THz sources, the mechanism of THz 
emission raises basic physics issues. To date most theory 
and experiments on coherent THz emission in air have 
considered localized plasmas with volumes ~ ^m'^ 
much less than Aipjjj. for typical THz wavelengths. In 
contrast, light string induced plasmas can be significantly 
longer than Athz along the propagation direction of the 
exciting laser meaning that coherent emissions from dif- 
ferent locations in the plasma need not add constructively 
[lo| . In this Letter, we propose incoherent light emission 
by a two-component plasma as a viable mechanism for 
THz emission since this does not rely on the establish- 
ment of relative phases between emissions from different 
locations in the plasma. We develop a microscopic the- 
ory which includes the Coulomb interaction between the 
charge carriers and the light-matter interaction with a 
quantized light field. First we show that even a homo- 
geneous plasma can produce broad-band radiation in the 
form of luminescence due to the electron-ion scattering. 
Second we make a mode function analysis for the specific 



rod-like geometry for light string induced plasmas in air 
in order to explore the angular dependence of the THz 
emission. 

Our starting point is the total Hamiltonian 



kill 



H 



A-p ' 



(1) 



where A = e,i sums over electron and ion contri- 
butions, and we next explain the meaning of each 
of these terms. This form of the total Hamilto- 
nian arises from the minimal substitution Hamiltonian 
^min = / d3r4'l(r)iE^^i^^^'A(r) apphed to our sys- 
tem where 5'A(r) and 'l'^(r) are the annihilation 
and creation operators for electrons and ions, respec- 
tively, m\ and qx are the particle masses and charges, 
A(r) = Z]q.o-(^«^q,cre''* '"eq,(T + h.c.) is the quantized 
vector potential operator. Here Aq = yJh/{2eoU}q V), 
and we introduced the boson operators Bq^a-{Bl^^^) 
which annihilate (create) photons with wave number q 
and polarization direction eq^u, in terms of which the 
free-field Hamiltonian appearing in Eq. is 7?l = 
J2q a- (-^q,o-^q,o- + I)' ^6 work in the A • p-picture 
and formulate our equations in momentum space, the re- 
lation between real and momentum space operators be- 
ing 5* a (r) = l/\/V'^^^e^^''^ax^i^, where the operator a^.k 
annihilates a particle of species A with momentum hk.. 
The kinetic energy term in the Hamiltonian is then 
given by H^^^ 

become 



k'^-^.kj with free-particle ener- 
A _ hj^^ g^j^j^ ^j^g light-matter interaction terms 



^A P - - X! "^k • Aq a^_^q/2"^^k-q/2, 
k,q 



(2) 
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with the canonical current matrix element = -^Kk, 
and the Fourier transformed operator of the vector po- 
tential Aq, along with the nonlinear ^^-contribution 



H^2 = 



k,q,q ,s 



Aq' • Aqa^ k+q,s«A,k+q%s. (3) 



Finally, Coulomb interactions between charges are incor- 
porated in the total Hamiltonian |^ using 



He 



= 1 Y v^-^' 
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^l,k«A'.k'«A',k'+qaA,k-q- (4) 



A, A'' 
k,k',q 



with the Coulomb matrix element V^^^'^ 

We work in the Heisenberg picture where the quan- 
tum average of an operator O evolves according to 
ihd{0) / dt ~ {[O^Htot]). Here, the many-body Coulomb 
and the light-matter interaction introduce a hierarchy 
problem such that e.g. photon numbers {B'^ B) are cou- 
pled to mixed correlations of {B'^a^a) which in turn 
become coupled to higher-order expectation values like 
{B'^a^a^aa). We truncate this hierarchy using a cluster 
expansion [l^ where multi-particle or mixed expec- 
tation values are consistently factorized. For example, 
photon correlation functions are split according to 

K',.'Sq,^) = K',.')(Sq..) + A(i?;,_,,Bq,.), (5) 

with (Bq.cr) related to the coherent component of the field 
and A(i3q, ^iB^^^r) describing quantum fluctuations. 

Here we are interested in incoherent light emission un- 
der quasi-stationary conditions after the exciting laser 
has passed, and an appropriate measure of the photo- 
luminescence is the total rate of energy emitted into a 
single mode (q, a) given by 



dt' 



Bl^B^^a) = r^,^A(i3t_,i?q^.), (6) 



where we have neglected any coherent field contribution. 
The Heisenberg equation of motion for the photon num- 
ber is found to be 



^A(St,i?q,,) 



q,o-'2k-q/2'^k-|-q/2; 



L k 



(7) 



where we have defined ek = cjA=e,k for notational sim- 
plicity. Here we have neglected the ^^-contribution 
in order to investigate the general emission mecha- 
nism. We will show later that its influence can be ac- 
counted for by generalized mode functions of the quan- 
tized field such that all our calculations remain valid. 
From Eq. Q we see that the incoherent emission of 
photons is driven by photon-assisted densities of the 



form A(_B^^ej^_q^2Sk+q/2)i which correspond to pro- 
cesses where an electron changes its momentum, while 
emitting a photon with the corresponding momentum 
difference. Under the assumption of large ion masses 
we have neglected a similar contribution from photon- 
assisted ion densities A{B^p^p). To evaluate the dynam- 
ics of A(i?^ (T^k-q/2^k+q/2) wc apply the long wavelength 
and low density limit; the procedure leads to 

«^7)7^(^q,^4-q/2ek+q/2> = 



2. rt'.k' (9) A«.<'-q/2ek'+q/2> 

k' 

'■Aq Cq.CT ■ Jk/kJ 



(8) 



where the last term corresponds to the spontaneous emis- 
sion source term. 

It is well known that a single-component, non- 
interacting plasma cannot emit or absorb radiation while 
simultaneously conserving energy and momentum. How- 
ever, in Coulomb interacting systems such restrictions 
do not apply. In the following, we concentrate on the 
scattering between electrons and ions, symbolized by F^ 
and F*^^ in Eq. ||HJ). We numerically confirmed that 
emission due to electron-electron scattering is negligi- 
ble compared to the electron-ion scattering. Additional 
mechanisms such as electron-neutral scattering will in- 
crease the scattering rates but are not expected to change 
the qualitative picture. A microscopic description of the 
scattering is obtained by including the coupling of the 
photon-assisted densities to photon-assisted correlations 
of the form A(i3^e^p^pe), deriving the Heisenberg equa- 
tion of motion for this quantity, and factorizing the re- 
sulting higher-order correlations in products of the form 
/(I — /)(! — /)A(i?^e'l'e). By Fourier transformation this 
equation can be formally solved and inserted into Eq. ((HJ 
in the form of frequency dependent scattering matrices 

(^*)^|k-k'| Ek" /k" r -,* , , 

^i^Mil)^— — + i^i ^ -^gJ ' (9) 



iS 



with diagonal elements F^(q) — X^k' ^k'^('?)' where 
W^|k-k'| denotes the statically screened Coulomb matrix 
element and large ion mass and low density limit 
are applied. The frequency dependence of the scatter- 
ing derives from the presence of the photon operator in 
the photon assisted correlations. Inspection of the de- 
nominator in Eq. shows that an electron changes mo- 
mentum due to Coulomb interaction with an ion. The 
corresponding change in kinetic energy is accounted for 
by photon emission while the total change of momentum 
of the electron is absorbed by the heavy ions. 

Figure ^ shows the calculated /q ^ as function of fre- 
quency for a homogeneous plasma with different elec- 
tron densities and temperatures. The temperatures lie in 



3 



40 



30 



10 



's 








.A 








r\ 
. ■ \ 
I- ^ 


Q. 

> 


\ 

\ 




. 1-. \ 
1 ■ \ 




\ 




\ ■ ^ 

' \ ■ ^ 
\ "■ ^ 






■x 

■ — 


\\\ \ 




1 


1 1 


100 


200 300 


- \\ ^ 




to 


[rad/ps] 










■ / '^^v 

i' ■ 1 









0) [rad/ps] 



FIG. 1: Luminescence signal Jq.^ 
plasma of density n 
(dotted), and n = 4 x 10 



^« c 

-|16 . 



as function of ujq — cq for a 
^ (solid), n = 2 X lO^'^ cm"^ 
(dashed) at T = 3000 K, 



and for n = 10^" cm'-^' at T = 2000 K (dash-dotted). The 
shaded area is obtained by multiplication of the solid line 
with the mode strength |uqp. Inset: data for n — 10^^ cm~^ 
and T = 3000 K after multiplication with uj^ with (solid) and 
without (dashed) frequency dependent scattering. 



the range of typical excess energies of about 0.25 eV for 
multi-photon ionization. We observe that the spectra are 
sharply peaked around uj — 0. The emission at a; = is 
mainly determined by the plasma temperature while the 
width of the emission tail increases with elevated density. 
In order to understand these results, we approximate the 
microscopic (fc-dependent) scattering matrices by an ef- 
fective constant damping rate 7off. In that case, Eq. (jH)) 
can be formally solved and inserted into Eq. 0, which 
upon assuming Maxwell-Boltzmann distributions for the 
electrons and taking the low density limit yields 



/PL ^ 



-'pi 



7off 



(10) 



If we take the main feature of Eq. (|5J) and assume that 
7cff increases linearly with density, we see that the high 
frequency luminescence is proportional to uj'^^jeS 
while the low frequency behaviour is determined by 
^pi/7cff 1 and thus independent of density. For low fre- 
quencies, Fig.^also exhibits the trend of decreasing emis- 
sion strength with decreasing temperature. The observed 
insensitivity of the high-frequency spectrum on temper- 
ature, however, cannot be explained with Eq. H10|l . Nev- 
ertheless, this equation can be used to fit an estimate for 
typical scattering times. For the case of n = 10^^ cm^"^ 
and T — 3000 K we obtain a scattering time of the order 
of 2.5 ps. 

Another shortcoming of Eq. H10|) lies in the fact that for 
high frequencies it is proportional to l/w^ such that the 
integrated power density proportional to ^ becomes 
divergent. In contrast, the microscopic, frequency depen- 
dent scattering leads to an exponential decay of 
as shown in the inset of Fig. ^ In order to obtain an 
estimate for the conversion efficiency of emission into the 



THz regime, we integrate the total emitted power density 
w2(27r2c3)-i/PL from w = up to w = 50 THz and ob- 
tain a value of roughly Pthz ~ 5x 10"* W/nr^. Comparing 
the total energy emitted by a plasma channel of radius 
i? = 30 /im and length of Z = 1 m over the characteristic 
lifetime t ^ Ins, E = TrR^lrVTaT. ~ 1-4 x lO^^^ J, with 
typical light-string excitation energies in the range of 
10 mJ, we obtain a conversion efficiency of 77 « 1.4xl0~^^ 
comparable to estimates for coherent emission mecha- 
nisms 

So far, we have completely neglected the influence of 
the ^^-contribution to the Hamiltonian. Instead of in- 
cluding the resulting terms into the equations of mo- 
tion, we pursue a different approach. We approximate 
the density operator in Eq. Q by its factorized expec- 
tation and assume that the overall density profile varies 
slowly enough to be taken as constant. In that case, 
is quadratic in the field in analogy to H]^ . The com- 
bined i/^2 +//l can be diagonalized by choosing the light 
modes according to the generalized Helmholtz equation 

V X V X Uq^<,(r) + Lol^{r) Uq,^(r) = uj"^ Uq,^(r) (11) 

with the space-dependent plasma frequency Wpi(r) — 

'^sam'^ • "^^^ computation of the new mode functions 
is a purely classical problem, the quantum character of 
the light field being fully contained in the photon an- 
nihilation and creation operators. If the extension of 
the plasma is much larger compared to characteristic 
electronic length scales, the charge carriers still see an 
approximately three-dimensional environment. In that 
case, we find only one modification that the lumines- 
cence spectrum is multiplied by a frequency-dependent 
mode function. For a thin slab of constant density with 
a thickness d well below typical emission wavelengths, 
this prefactor for emission perpendicular to the plane of 
the plasma is simply given by 



|u,(r = 0)|- 



I 4 g2c4 



Thus, we see how frequencies well below the plasma fre- 
quency are effectively blocked out. The shaded area in 
Fig.^shows the resulting spectrum for an assumed thick- 
ness d = 3 /im. For thicknesses d 3> c/ujp\ all frequencies 
below the plasma frequencies would be completely absent 
in the emission. 

The geometry of a laser induced light string enters 
our computations in exactly the same way. The eigen- 
modes have to be computed according to Eq. ((TT|l in 
the presence of the spatial density profile 71 (r) of the 
plasma column^^. Our problem has well separated 
length scales, namely, the thermal wavelength of the elec- 
trons Ath — 1 nm is much smaller than the diameter of 
the plasma column i? « 30 /im, which is in turn much 
smaller than typical THz wavelengths Athz- As far as 
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FIG. 2: Mode strength \uq^a{r = 0)^ for (a) cq = lps~^ < 
LOpi and (b) cq ~ 4.5 ps~^ = O.Suipi. The corresponding mode 
strength in vacuum would be 1. 

optical coupling is concerned, the exact form of n(r) is 
therefore relatively unimportant as long as it is strongly 
localized, say at the origin of the {x — y) plane. We have 
solved for the eigenmodes of Eq. by using a trans- 
fer matrix technique, extending the approach of Ref. 
to oblique incidence. All solutions can be classified by 
the wave vector q of an incoming plane wave and its po- 
larization a distinguishing between TM-mode (incident 
electric field in the plane of incidence) and TE-mode (in- 
cident electric field perpendicular to electron string). 

The mode strength |uq,j(r = 0)p for the TM-mode 
is presented in Fig. |21 for the case of a step-like density 
profile with a carrier density of n = 10^^ cm~'^ over a 
cylinder with radius i? = 30 /im. For frequencies well 
below the plasma frequency, we find a dipole radiation 
pattern qualitatively similar to that observed experimen- 
tally in Ref. Q. The main emission for those frequen- 
cies is perpendicular to the plasma rod aligned along the 
horizontal axis of Fig. |2] This emission pattern can be 
traced back to the increasing suppression of light modes 
below the plasma frequency as the emission angle de- 
viates more from the normal direction. For frequencies 
around and above the plasma frequency, the angle de- 
pendence drastically changes since this suppression does 
not exist anymore. For certain frequencies and emission 
angles resonance effects can be observed which depend on 
the precise choice of parameters. In the given example, 
we observe a preferred emission angle around 40 degrees. 
The mode strength of the TE-mode does not strongly 
depend on the angle of incidence since for any angle the 
electric field component lies perpendicular to the light 
string. 

In summary, we presented a microscopic theory for the 
photoluminescence emitted from fs laser-induced plas- 
mas, and have shown that even a homogeneous two- 
component plasma can emit due to electron-ion scatter- 
ing. The angle dependence of the THz emission for light 
string induced plasmas was explored using modal anal- 
ysis, a dipole pattern for emission frequencies below the 
plasma frequency being found in qualitative agreement 



with the recent experiment |3. In contrast, a drastically 
different angle dependence for emission frequencies above 
the plasma frequency is predicted. Finally, we note that 
although our analysis was motivated by studies of THz 
emission in air |^ |^ Q , our results apply more gener- 
ally, for example, to glasses or electron- hole plasmas in 
semiconductors . 
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